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ABSTRACT 

Transmission of the power is the major challenge as utilities located at the remote places, hence proper 
functioning of power transformer is to be maintained. Many times, transformer encounters with no-fault condition such 
as inrush current, which results in unwanted tripping of transformer. So, it is necessary to discriminate fault and no¬ 
fault condition on transformer, so that unwanted tripping of circuit breaker can be avoided and uninterrupted service 
can be maintained. 

This paper, present a technique based on the Continuous Wavelet Transform (CWT) to differentiate among 
faulted and un-faulted condition on power transformer. Transformer modeling and simulation of faulty and healthy 
condition are carried out using MATLAB/Simulink software. By employing CWT, suitable features are extracted as the 
waveform resulted from the faulty condition, and healthy condition are different in nature, after first zero crossing of 
fault current and Inrush current. 
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INTRODUCTION 

In recent years, many researchers in mathematics and signal processing have developed wavelet 
transform for analyzing signals in time and frequency domain. This transform is most suitable for non- stationary 
signals. Robert et al (1996) proposed the use of wavelet transform for analyzing electromagnetic transients 
associated with power system faults and switching. The frequency domain data produced by wavelet transform is 
useful for analyzing source of transients. The choice of mother wavelet for different power system applications 
have been proposed by santoso et al (1996). According to him, for short and fast transients, db2, db4, etc. are 
preferred and for slow transients, db8, dblO shows good results. Use of wavelet for Time resolution for high 
frequency components of fault transients is discussed by Fernando and Ali Abur (1998). The methods for 
calculating effective values of voltage and current in wavelet have been worked out by Zhu (2004). 

Many techniques have been proposed to identify internal fault current and inrush current, these 
techniques were based on signal processing using fourier transform method, calculation of time interval between 
two peaks of fault current and two peaks of inrush current, calculation of peaks for both currents and comparing 
them to develop an appropriate algorithm or model. In this paper, 3-phase ground fault has been considered as 
faulted condition and Inrush current is taken as no fault condition on transformer [10]. 
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Under abnormal condition, power transformers are protected by differential relay. This relays are often get tripped 
due to inrush current, which is developed during charging of the transformers. Hence, treatment of inrush current as fault 
current is the main drawback of the protection system. So, there is need to develop a system, which can differentiate 
between faulty i.e. internal faults and no fault i.e. inrush current. If the sensing of inrush current as fault current by relay 
will make the circuit breaker to operate, which is unreliable, since a large area is disconnected from being supplied. The 
malfunctioning of relay results in unreliable power system. Hence, the property of inrush current and fault current can be 
used to discriminate between them and stability of the power system is improved. 

The mathematical tool called, wavelet is becoming popular method for condition monitoring of large 
transformers. Abnormal condition on power transformer has to be checked by disconnecting it, and then opening each part 
to locate the fault. This traditional method is time consuming and very expensive, also disconnecting it from utility is not 
good choice. 

Hence, Wavelet Transform is used to analyze the signals, so that exact location of fault can be identified without 
opening or disconnecting power transformer and online condition monitoring can be done. 

FAULTS ON TRANSFORMER 

Although transformer is static device, due to abnormal system conditions, internal stresses arises. Generally, 
transformer is encountered with faults mentioned below 

• Over current due to increased loads and external short circuits 

• Fault at the terminal 

• Faults in the winding 

• Incipient faults 

Thermal and mechanical stresses results from these faults in the winding of transformer and its terminals. 
Overheating is the result of thermal stresses, which ultimately affects the transformer by affecting its insulation system. 
Insulation deterioration leads to the development of faults in the winding. Sometimes, the cooling system of the 
transformer fails and results in overheating, hence safety of transformer is very important. The reactance of transformer 
helps in limiting the short circuit current of the electrical transformer and short circuit current is high at low reactance. The 
generally occurring winding faults are earth fault and inter turn fault. The incipient faults are nothing but the faults, which 
have no instant effect but, the persistence of this fault for long time may lead to the major fault. For example, the incipient 
faults can be due to damage of insulation between core and lamination, reduction in level of oil due to leakage, stoppage of 
oil flow paths. As a result, transformer overheats which may damage the insulation, hence resulting into short circuit [1, 6]. 

INRUSH CURRENT 

It is an instantaneous current, drawn by transformer, when operating under no-load during starting. When 
transformer is initially energized, a large current 10 to 15 times the rated transformer current flows for several cycles. The 
inrush current is the result of hysteresis property of material as de-energizing a transformer retains some magnetism, which 
develops residual flux and when the transformer energizes, this residual flux voltage helps in developing a transient current 
of large magnitude known as Inrush Current. This current is not the fault current, as it flows only for few cycles and does 
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• capacitance of winding 

• characteristics of transformer core material 


• current chopping characteristics of circuit breaker 

• core gap factor etc. 

The magnitude and the time length of inrush current depends on transformer size, power system size, switching 
instant, material of core of transformer [13]. 

The equation for generation of inrush current is given in (1), where a three phase transformer is connected to AC 
supply. V and A are the instantaneous flux in the core and voltage drop across the primary winding. 

dA 

y = — a) 

dt 


We can say from equation (1) that flux wave form is the integral of voltage wave form and in operating 

transformer these are shifted 90 . For energizing transformer, the flux initially starts from zero point, so voltage starts 
from zero point, flux will reach approximately twice its peak as it integrates the area under the voltage curve, for first half 
cycle, leading to the saturation of transformer. During saturation, disproportionate amount of flux has to be provided, 
which is meet by mmf produced by winding current in the core, and hence flux becomes double of its normal peak. Thus, 
generates nearly ten times large current know as magnetizing inrush current. 


The general equation of magnetizing inrush current can be expressed as 


i(t) = 


Z, 


*K W *K^ *(sin(ox-</>)-e 1 .sinor) 


( 2 ) 


Where, 


V m : Maximum applied voltage; 

Z t : Total impedance under inrush, including system; 
(/>: Energization angle; 


t: Time; t 0 : Point at which core saturates; 

T: Time constant of transformer winding under inrush conditions; 
CC : function of 1 0 ; 


K w : Accounts for 3 phase winding connection; 


K t : Accounts for short-circuit power of network. 
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Figure 1. shows inrush current resulted from simulation in MATLAB for the proposed system for 1ms duration. 

SYSTEM DISCRETION 

The system under study is composed of a 3 phase synchronous generator delivering power to load through 3 phase 
power transformer. A 3 phase transformer protected differential relay is by employing two circuit breakers on primary as 
well as secondary. 

The difference between the currents of two CBs, which are the outputs of primary as well as secondary is given to 
the proposed relay. The developed system is composed of 

• Three phase 3000MVA power source, 120kV 

• Two- three phase circuit breakers on the transformer sides. 

• Three phase star/star 47MVA, 120/25 KV power transformer. 

• 3-30 KM feeder. 

• Three phase load. 

• The required current transformers 

As it is known that fault current is cause of residual magnetism in the core, hence core parameters were adjusted 
to obtain inrush current, and fault current can be obtained by three phase ground fault in matlab/Simulink model. 

Table 1 shows transformer parameters used for the simulation. 


Table 1: Simulated Transformer Main Parameters 


Parameters 

Ratings 

Transformer Connection 

Y/Y 

Rated power 

47MVA 

Voltage Ratio 

120/25kV 

Rated Frequency 

50 Hz 

Primary winding Resistance/phase 

0.003 pu 

Primary winding Inductance/phase 

0.09 pu 

Secondary winding Resistance/phase 

0.003 pu 

Secondary winding Inductance/phase 

0.09 pu 

Core resistance 

300 pu 

Magnetizing Reactance 

300 pu 


In this model, Wavelet based discrimination technique is used to distinguish between fault and no fault situation 
in transformer. 

SIMULATION OF FAULT CURRENT AND INRUSH CURRENT 

The developed model is used to carry out simulation of fault current as well as inrush current. Since inrush current 
is result of residual magnetism, hence in modeling of transformer, its core saturation behavior is considered. Simulations 
results for the two currents are shown in Figure 2.and Figure 3. From the result, we can see that both waveforms are 
different in nature, in case of inrush current, the magnitude is increasing abruptly for few seconds and decreasing at fast 
rate, and in case of fault current, magnitude increases for few cycles, but after its magnitude remains constant if not 
interrupted by circuit breaker. Developed star/star transformer matlab/Simulink of the power system model is shown in 
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Figure 1: Developed Star/Star Transformer Matlab/Simulink of Power System Model 

The discrimination between fault current and inrush current can be done by using the characteristics of both 
curves, as inrush current curve has small slope at starting and as it progress its slope, gradually increases. Moreover, DC 
components in inrush current are more and last for few cycles, but in case of fault current, the slope of the curve is more at 
every situation and DC components are very rare. Though the difference in both cases is very small, but they are sufficient 
to distinguish between fault and no fault condition of transformer. 

The early method of discrimination i.e. differential protection was based on harmonic content of differential 
current. In which, the ratio of second harmonic to the fundamental component of differential current for fault current would 
be less than inrush current. Hence, clearing of faulty condition has been successful, but the main drawback of this method 
was the harmonic production, which is also possible due to saturation of CTs and long transmission lines, so, there was 
possibility of operating differential relay under no fault condition, also hence this method is unreliable. 

The other methods used were based on fluctuating waveforms of differential relay and time difference between 
two peaks of internal fault current and inrush current. From studies, it has been proved that difference in time for fault 
current is less than inrush current. Wavelet based discrimination technique has shown great advantage in identifying as 
well as clearing fault when occurred, as it is very effective and fast technique which can clearly recognize difference 
between fault and no-fault condition on the transformer. In this technique, different samples can be obtained for few cycles 
of fault as well as inrush current. Then, this data can be compared in terms of their magnitude and then can be used to 
actuate the differential relay. 
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Figure 2: First Four Cycles of Internal Fault Current 



Figure 3: First Four Cycles of Inrush Current 


WAVELET TRANSFORM 


A mathematical tool for analyzing the signal is Wavelet Transform. It is the extension of Fourier Transform 
method. This method provides the information regarding the frequency component present in the signal at particular time. 

There are two types of wavelet transform, Continuous Wavelet Transform (CWT) and Discrete Wavelet 
Transform (DWT). CWT uses inner products to measure the similarity between a signal and an analyzing function. It is 
defined as, the sum over all time of the signal multiplied by scaled, shifted versions of the wavelet function. Calculating 
wavelet coefficients at every possible scale is a fair amount of work, and it generates an awful lot of data. In DWT, we 
choose only a subset of scales and positions at which to make our calculations. 

A mathematical tool for analyzing the signal is Wavelet Transform. It is the extension of Fourier Transform 
Method. This method provides the information regarding the frequency component present in the signal at particular time. 

There are two types of wavelet transform, Continuous Wavelet Transform (CWT) and Discrete Wavelet 
Transform (DWT). CWT uses inner products to measure the similarity between a signal and an analyzing function. It is 
defined as the sum over all time of the signal, multiplied by scaled, shifted versions of the wavelet function. Calculating 
wavelet coefficients at every possible scale is a fair amount of work, and it generates an awful lot of data. In DWT, we 
choose only a subset of scales and positions at which to make our calculations. 

CONTINUOUS WAVELET TRANSFORM 


Wavelet means a small wave function and compactly supported. The standard function for generating the other 
window functions is mother wavelet. At lower frequencies WT has excellent frequency and poorer time resolution, and at 
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higher frequencies the WT has excellent time and poorer frequency resolution. 

The continuous signal x(t) can be expressed in Wavelet Transform as 

wmr) = -M*(0*g(— )dt (3) 

A — T 

Where, g( -) is an analogous function used to transform signal x(t). The characteristics of analyzing function 

a 

g (t) is that it must be short and oscillatory. 

CC- Dilation Factor or scale, in terms of frequency, global details of a signal can be acquired from low 

t — T 

frequencies and information regarding hidden pattern can be obtained from high frequencies. Mother wavelet g( -) 

a 

can be expanded or contracted in time depending on whether 0C> \ or Cfcl.If CC> 1 then signal is expanded in 
time and if CC < 1 then it is contracted in time. 

t = translational factor, gives information of time. 

Various mother wavelets have been derived such as Daubechies wavelet, Haar wavelet, Symlet, Biorthogonal, 
Coiflet for analysis of signal into a set of approximate and detailed coefficients. It uses MRA (Multi Resolution Analysis) 
method to extract the information of the signal [4, 10]. 

DAUBECHIES WAVELET 

The explicit function of dblis Haar Wavelet other than these Daubechies (db) wavelet has no exact expression. 
The db can also be written as dbN where N is the order, db having Supporting length of if/ and (j) is 2N-1. The number of 
vanishing moments are more, hence analysis of signal is made more clear. 



Figure 4: Daubechies Wavelet 

N is also the number of vanishing moments in y/ . db wavelets are not symmetrical and most of the functions are 
regular at certain points. 
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SIMULATION RESULTS AFTER APPLYING CWT 


The proposed system in Figure 1 is simulated in MATLAB/Simulink and results are shown in Figure 5 and Figure 
6.for daubechies wavelet (db4), applied at level 5 for the fault current and inrush current. The reason of choosing db4 is 
that the fault and inrush current are short and fast signals, and according to santoso et al (1996), for such signals, db2 and 
db4 are preferred so as to extract more information. 

Different levels of wavelets applied on the internal fault current obtained from simulation are shown in Figure 5. 
The frequencies generated due to fault can be seen. The 5 levels show different characteristics and coefficients are also 
different. The time span taken for internal fault is t =11 ms and resultant frequency (db5-dbl) is shown in Figure 5 and 
Figure 6. For comparison purpose, we chose dbl coefficients and from the Figure 7, the characteristics are compared. In 
the same way, the fault current is restricted to only for four cycles, and the analysis is made and frequency levels from 
wavelet are demonstrated. 


Wavelet coefficients for Fault Current 
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Figure 5: Wavelet Coefficients for Internal Fault Current 
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Wavelet coefficients shown in Figure 6.are also obtained by applying different levels of daubechies wavelet for 
analyzing the inrush current resulted from simulation of the proposed model as shown in Figure 1. 

Wavelet coefficients for Inrush Current 



Figure 6: Wavelet Coefficients for Inrush Current 

The wavelet coefficients shown in Figure 5.and Figure 6.are arranged as db5, db4, db3, db2, dbl in descending 

order. 

Figure 7 presents the profile of db4 applied for level 1 for internal fault current and magnetizing inrush current 
respectively, and is used for discrimination purpose. Few values of coefficients can be taken for comparison of both the 
conditions. 60 coefficients are taken for discrimination purpose and the range can be seen in Figure 7. dl components are 
considered for discrimination purpose, so that a clear difference in magnitude and frequency is obtained. The difference in 
both the signals is observed only for first four cycles, since the magnitude of inrush current is more for first few cycles. 
Different profiles can be obtained by considering d5, d4, d3, d2, and dl. This technique can be validated considering the 
signatures of both internal fault current and inrush current after first zero crossing. 
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ell for Fault Current 




Figure 7: D1 WT Coefficients for Fault Current and Inrush Current 


The analysis is made for values of dbl wavelet coefficients and these coefficients for internal fault current are 
more than the corresponding values for inrush current and this is given for one half cycle in Figure 8.and Figure 9. 

Also, the slope of inrush current is very less initially, and gradually increasing reaching the next cycle, but the 
slope of fault current is sharper. 
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Figure 8: Internal Fault Current after First Zero Crossing 



Figure 9: Inrush Current after First Zero Crossing 

So, after first zero crossing as shown in Figure 8.and Figure 9, it can be observed that values obtained from both 
wavelet coefficients varies drastically. After zero crossing, the values of wavelet coefficients in case of inrush current are 
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first negative, then after quarter cycle they are positive, but for internal fault current they are always negative. These values 
are used to discriminate between internal fault current and inrush current in the power transformer. These values are 
sufficient to make the relay sense, the fault correctly by making decision that, inrush currents are not the faults. The values 
obtained from this technique can be used as input to the digital relay for decision making purpose, so that unwanted 
tripping is avoided and uninterrupted supply can be provided during no-fault condition. This will not only make the power 
system efficient, but life of protection system also get enhanced. For implementation purpose, the initial cost of the power 
system increases, but running cost will drastically come down because of proper functioning of differential relay. 


dl for Fault Current 



Figure 10: Elements of dl Wavelet Coefficients for Fault Current 


dl for Inrush Current 



Figure 11: Elements of dl Wavelet Coefficients for Inrush Current 

Figure 10 and Figure 11 shows the elements of dl wavelet coefficients obtained for both cases and are sufficiently 
different to discriminate internal fault current and magnetizing inrush current, and these values can be given to digital 
relays. 

CONCLUSIONS 

We studied different types of faults occurring on transformer and characteristics of inrush current. The brief study 
of wavelet transform is made along with daubechies wavelet, so that analysis of the signal can be understood clearly. Also, 
the simulation of internal fault current and magnetizing inrush current for developed model of power transformer has been 
done and coefficients are extracted. Differences between the magnitudes of wavelet transform coefficients in specific band 
of frequency, obtained from internal fault and inrush currents. This element of dl wavelet coefficients are used for 
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discrimination purpose, as they are sufficiently different in both cases of fault and no-fault condition. The wavelet 
transform method has been proved to be the efficient method and fast technique to discriminate transformer internal fault 
and no fault condition. This can be again simulated in Fuzzy Logic. 
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